Abstract 


This registration report details the design, 
construction and software support for a digital torsion 
meter capable of recording transient twist of any shaft, 
Graphical results are presented for the delta start 
transient for an induction motor driven electromechanical 
rig and compared to results generated by a lumped 
parameter model for the system. Indications of other 
tests that will be carried out and improvements to the 
software and hardware are given, It is projected that the 
accuracy of this device may be improved to 0.02 degrees 
and the signal reveal a torsional signature of the states 
of operation and health of a machine, 
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Graph 1: Velocity - raw data 


No error detection and the minimum of calculation to 
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The routines for the detection and correction are 
applied to the data for graph 1 


Graph 3: Velocity - as previous and no over-twist 


The routines for the detection and correction of over 
twist are applied to the data 
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The routine for generating the scaling factors is 
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subsequent processing 
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Symbols 
angular acceleration 
compliance of the belt transmission 
coefficient of the frictional torque acting on the 
rotor 
d/dt 
Young’s modulus for the belt 


function relating angular displacement to the time 
after start 


function for the rt state variable of lumped 
parameter model 


tension in the belt arms 


Suffices 
t taut arm 
s slack arm 


the inverse function of f[{t} 


rotational inductance of the generator 
intermediate results used for Runge-Kutta algorithm 
Suffices 


ryn {=<r=<4, 1=<n=<number of variables 


iteration interval of the Runge-Kutta algorithm 


current in an equivalent circuit of the motor = or 


generator 
Suffices 
hts ee full model of motor 
K1,y4 two coil model - stator 
x2,y2 ey - rotor 
0,Q@ D-Q model stator 
d,q a rotor 
current vector for the motor 
Suffices 
3ph three phase vector 
2ph two phase vector 
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moment of inertia 


Suffices 
4 motor rotor, coupling and half shaft 
2 coupling and pulley 


S generator and pulley 
currents in generator 


Suffices 
*, armature coils 
= field coils 


integer count of the number of transitions 


scaling factors for the compensation of strips 


Suffices 

4 for the velocity measurement 

2 for the torsion measurement 
stiffness 

SH shaft 

b belt 
mutual inductance between coils of the motor 

Suffices 

u,v O<u, V7, UddV 

full motor model 
m D-Q model, mean level 


self inductance of coils 
Suffices 
Aichteee full motor mode! 


number of transitions per revolution at 
measurements are taken 


resistance of coils 
Suffices 
Pree rcorl) full motor model 


resistance of various components 


Suffices 
L generator load 
R rotor bars 


Ss stator coils 
radius of pulleys 


angular displacement 


Suffices 
1 general 
2 


MS motor end of shaft 
sp pulley end of shaft 
G generator 
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24 time elapsed 
Suffices 
1 interval of encoder period 
2 time between transitions at shaft ends 
iil torque 
Suffices 
B transmitted by belt 
diff }0SS of torque across the belts 
£4 frictional torque on motor 
#2 frictional torque on generator 
G generator load torque 
M motor torque (airgap) 


PG torque on the generator due to belt 
pre initial tension of belt 


v phase voltages across motor 
Suffices 
1,...,6 full model 
angular velocity 


xp rth state variable of lumped parameter model 


Ww mains supply angular frequency 
Lh coefficient of friction between belt and pulley 
p angle of contact between belt and pulley 
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Chapter 1 
Introduction 
4 introductory 


"The need to measure accurately the torque, 


or torsion moment, in ae shaft, ... has 
exercised the ingenuit of engineers for 
many years.” 


Ford and Douglas wrote this to introduce their paper 
in 1949, aware of the fifty years of development that had 
preceded. In the forty years since yet more sensitive and 
accurate meters applying the currently available techno- 
logies have been built. This research is not to develop 
yet another torque meter, but to monitor -the twist by a 
method that will allow application to a number of 
rotating structures and generate signals to indicate the 


states of operation and health. 


The research began as a study into the feasibility 
of condition monitoring a erankshaft while in service, 
it was taken that the feasibility would be proven by the 
development of such a technique. This technique would 
measure the torsion induced in the shaft due to the 
torque transmitted and be able to pick out transients 


and characteristic twisting vibrations, 


‘Ford, Hugh and Douglas, Alan; "The measurement of 
torque in shafts"; Engineering; 27/5/1949; pp 481 - 
484, 505 - 507 
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There are many uses for such a_ technique, however, 
the initial stages are being carried out on ane electro- 
mechanical rig. This: is to verify the developed 


technique against previous results and _ theory. 


1.2 Review 

Many methods exist to measure the torque and power 
output of a prime mover. For example, some monitor a 
standard section of the transmission system®. There are 
other systems developed more recentiy that will monitor 


the torsional vibrations of a rotor by laser 


techniques3, 


The former usually require access for attachments” to 
the transmission chain or even a standard section to be 
inserted, changing the characteristics of the measured 
system, The latter technique can only monitor the 


velocity at one point from which the torsional vibration 


2ingham, A. J.; "Torque measuring systems. Measurement 
Measurement from rotating shafts, associated 
electronics and calibration methods"; Industrial 


Sensors seminar; Advanced manufacturing techniques, 
demonstration heid at Regional electronics centre, 
University of Hull; 23/6/88 


3£astwood, P, G. and Halliwell, WN. A.; "Laser-based 
measurement of torsional vibration"; SPIE Optics in 
Engineering Management 599; 1985; pp 358 - 366 

ibid, "Laser Technology for Condition Monitoring"; 
COMADEM 88; Eds. Raj Rao and Anthony Hope; Conference 
at Birmingham Polytechnic; 19-21/9/88; pp 154 - 159 
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has to be deduced, 


Recent work from the Hitachi industrial laboratories 
measure the torque and power output by the mathematical 
inference from measurement of speed and the theoretical 
knowledge of the particular prime movers4, This technique 
- state observer - has the limitation that there must be 
considerable a priori Knowledge of the prime mover = and 


the system of which it is a part. 


1.3. Overall aims 

It is intended that the monitoring technique 
described below can give information about the particular 
system that is as yet not possible to measure, and this 
measurement may be carried out in a way that does not 
affect the characteristics of the rig. Also that the 
observation of the system is independent of Knowledge of 
the system, like the signals obtained from an 


accelerometer. 


The first stage of the project, the stage reported 


4onmae, Tsutomu; Matsuda, Toshihiko; Kanno, Minoru; 
Saito, Keiji and Sukegawa, Takashi; "A Microprocessor- 
Based Motor Speed regulator Using Fast-Response State 
Observer for Reduction of Torsion Vibration"; IEEE 
Trans. Industry Applications 1A-23 (S)% 
September /October 1987; pp 863 - 871 


Jti300988 4 


CNAA Registration report 


here, is to develop a monitoring technique. This includes 
analysis of the technique’s limitations and methods of 
deducing the state of the system from the measurement of 


torsion (and any prerequisite measurements). 


Also reported here is the development of the test 
rig mathematical model which was used in the first stages 


to verify the results of the measurement technique. 


The principle exploited and tested on this rig is 
that of described by Pratt®, The main difference being 
the use of the newer microcomputers to ‘carry out the 
processing and display. This report deals with the 
development of this on the drive shaft and the necessary 


computer software to support it. 


The ideas that form the basis of the research are, 
in summary:- 
A version of the angle of twist detector 


that is supported by a microcomputer to 
allow sophisticated analysis of its measure- 


ments, 

Spratt, JP; George W. and McMullin, Paul G.; "Opto- 
electronic apparatus to generate a pulse modulated 
signal indicative of the state of health of the 
mechanical state of a system”; United States Patent 
3871215; granted 18/3/1975. 

Pratt, Jr, George W.; "An Opto-Electronic Torquemeter 
for Engine Control"; Society of Automotive Engineers 


SAE paper 750070; 2/1976. 
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A dynamic model of the system that faith- 
fully follows torsion transients, 


1.4 Application to other systems 


When torque is transmitted through a rotating shaft, 
then the resulting twist along the shaft may be measured. 
By the nature of the opto-electronics used the meter 
is particularly non-invasive, that is it does not 
require alteration of the system measured other than 
suitable optical tracks to be attached to the trans- 
mission train and access made for the light beam. This 
system is useful in areas of limited space. By this, it 
is possible to use the system on any torque’ transmitting 


shaft. 
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CHAPTER 2 


The torsion meter principle of operation and execution 


2.4 Principle of operation 

When = any shaft is transmitting power, the torque 
twists that shaft. This meter measures this twist by 
determining the ratio between the time for two marks of 
Known angular displacement at one end of the shaft to 
pass a reference point and the time between the first 
mark mentioned above and a similar mark at the other end 
of the shaft. The first interval is inversely 
proportional _ to the angular velocity. The second interval 
gives a measure of the power if fn manke are initially 


coincident. Hence their ratio reveals the twist. 


If this operation is carried out frequently within 
a revolution and the times are measured accurately then 
the time history of the velocity, power and twist of the 
shaft may be displayed in detail. The mathematics 
supporting this method and its accuracy and sensitivity 


are worked out below. 


A note on terminology; the interval which reveals 
the velocity is referred to as the ‘period’, and that 


which reveals the power, the ‘interval’. 
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2.2 Mathematical Theory of Operation 


The period of the signal, ty, is used to measure 
the average angular velocity over that particular 


interval, 


(2.4) 
The displacement between the two code discs !s given by 
(82-84) =(ta) en, 
(t4)Ng 
(2.2) 
The above formulae are necessarily approximations, 
Usually it is taken that angular displacement is a 
function of time: 
@=f{t) 
(2.3) 


and 


De=Df{t}. 
(2.4) 


But for the purposes of this analysis it is more 


convenient to take the inverse function, g say, 


t=9(0} 
(2.6) 
and 
Gygltizg’(aj=t’. 
(2.6) 
This is to allow advantage to be taken of the fixed 
angle over which the measurements would be taken, But 
strictly the following analysis is only certain for 
unidirectional rotation, 
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The time for one period from one of the coding discs is, 


am(k+1) /ng 


t= t’de, 
1 
eTk/Ns 
(2.7) 
= glam(kt+1)/ng} - gletk/ng}. 
(2.8) 
This gives rise to the Taylor expansion, 
ty= 9g’ L2mk/ng} (21/ng) 
+1/a9’’ (amk/ng} (2m/ns) +... 
(2.9) 
.The first term gives the equation (2.1) with 
De = 47g? 


and and error due to the series truncation of the order 


of the series’ second term (if the series is convergent). 


Similarly for the twist, the time between successive 
leading edges is given by 
etk/ns eTk/Ns 


= g’ de, 


ie) te) 
(2.10) 


where h is the angle function of time for the other 
end of the shaft. Integrating out and using the following 
expansion 


to= hlewk/ng} - gl2tk/ng} + gfO) - htO}, 
(2.44) 
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defining 9(0),h(O0) as O, 
hl2tk/ng}=glewk/ng} +9’ (2mk/ng} (82-04) 
+1297? (2mK/ng} (82-04) 2+... 
(2.12) 
From (2.11) and (2.12) 


to= 9’ (2mK/ng} (82-84). 
(2.43) 


From which (22) follows and the error due to the 
approximation is of the order of the third term of 
(2,12). The value for g’ that is used in ae evaluation 
iS worked out from equation (2.1). This gives a 
cumulative error of 


\yatg?? Lemk/ng}} [(2m/ng) 4+ (82-04) 41172. 
(2.14) 


For most measurements the twist is less than a2t/ng then 
the total error has an upper bound of 


‘aig?’ Lemk/ng}! (2t/ng) @ 
(2.15) 


This requires an upper bound on the second 
derivative of time by angle to give an error independent 
of the functional relationship between time and angle. 
This derivative may be expressed in terms of angular 


velocity and acceleration. 


If v is the angular velocity and a is the angular 


acceleration then 


d*t = -a 
dee v3 


(2.16) 
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2.3 The construction of the discs 

The encoding discs used on this rig are of ABS 
plastic cut to a circumference of 480mm and mounted on 
each end coupling of the drive shaft&, To the outside 
edge of each disc is fixed a photographic strip of 


equally long black and white zones, 


Each zone is 12mm iong giving a total of twenty 
black/white transitions at equal angular intervals around 
the disc. To improve contrast the photographic film is 


mounted on an aluminium foil tape attached to the disc’, 


2.4 Flexibility in design 

The dimensions of the rig and early techniques 
developed to mark the discs constrained the discs to a 
circumference of 480 mm. This was to allow them to be 
constructed easily. However, there is little constraint 


on the size that they might be in a (final apparatus. 


The use of photographic techniques allows a = shorter 
strip with finer divisions up to the resolution of the 


interface system. 


Ssee Figure 1, Code dise location 
See Figure 2, Code disc construction 
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The signal processed by the computer is derived from 
the changes in reflectivity at the discs’ edges. 
Considering the reflectivity at one point above the edge 
of each of the discs gives a signal§ which does not 


follow the perfection of the above equations. 


The errors induced into the mathematical analysis by 
the real system stem from the imperfect nature of the 
production of the strips, The formulae assume that the 
transition occur at Precisely 2t/ng radian intervals 
about the circumference of the ANTS ‘This is not the 
case. The place of transition may in practice vary from 


its theoretical location by */_0,25mm on the surface. 


However, any errors in the measurement due to. this 
misalignment of the coding strips will be reflected in 
the measurements as a periodic variation over each 
revolution, This periodic variation may be compensated 
for by the use of scaling factors that are tuned to each 
zone length between transitions. This is achieved by 
running the discs at a constant speed and at a constant 
displacement and comparing the times for each transition 


with that calculated from the period of one revolution, 


8see Figure 3, Light signal 
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For the signal that is used for the velocity 
measurement 9: 
te) +g tee eee ST 
(2.47) 
where N is the number of transitions in one 
revolution and J 1s the period of that revolution 
So for each scaling factor ky!, 


kyity'=1/N 
(2.18) 


Similarly for the measurement of twist each period of 
the original signal giving the interval is given by 


to'tt-tolety!, 


(2.49) 
When the appropriate scaling factors are included, 
Kal tttoitt-Kaltolekyityt=T/N 
(2.20) 
From equation (2.18) this reduces to the ratio, 
tal*tyta! skal sea! +1 
Kase 1) 


The scaling factors of ali times, to, bear this ratio to 
each other. It is then defined that 


Ka=1 if ta=T/N 
(2.22) 


and all other factors are scaled to this. This reduces 
the fluctuations in the raw data before processing to a 
residual level which 1s) indistinguishable from = other 


sources of error. 


9see Figure 4, Signal Compensation 


JtI300988 20 


CNAA Registration report 


2.6 Driving the optoswitches and signal conditioning 


fis 1s a desemption of the Mark 1 interface 


circuit that drives the reflective opto-switches!? inter- 


faces the signals to the general interface port of the 
BBC microcomputer’s anctluding the Archimedes) usually 
termed the ‘User Fort’, A dagram of the eircurt iS 
displayed in Gircuit 1, The signals from the transistor 
follow the changes in the reflected light intensity from 
the -surface moving beneath it. The rest of the circuit 


compensates tor the ramped change in the signal with a 
second transistor and conditions it to TIL levels with a 


chan of Senmtt triggers. 


2.7 The microcomputer interface Mark 1, ’‘’User Port’ 


The information from the rig consists of an 
electrical signal. Ine computer has to turn ts signal 


into bytes of nrormation, 


The interface and communications from the rig were 
developed on a &BC microcomputer Master system and 
upgraded to an Acorn Archimedes. However, this first 


system employs the ‘User Port’ of these computers in the 


1ORS Reflective Optoswitch part no 307-913; Data sheet 
4276 3/3/87 
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same active way. The algorithm and use of the supporting 
component of ths first interface are discussed. 


Some of the draw backs in terms of the = lmitations 
and errors " data resulting from the use of this 


interface are outlined, 


A 6522 VIA supports the User PORT. The features that 


are particularly used by the conversion process are:- 


A) The generation of interrupt signals to the system 


CPU on active edges on signals apphed to its 
control lines and 
B) Two general purpose timers which may generate 


interrupts after a preset time interval 


To trace the signal through the computer system = and 
the software developed, let us consider Just one 
transition at the opto-switch, A leading edge appears on 


one of the control lines to the VIA. 


The internal registers of this chip would have been 


set to allow an interrupt to be sent to the CPU, 


The microprocessor detects the interrupt status, 


Jjeaves its current program and carmes out an interrupt 
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chosen on the basis of the VIA status. 


The routine for clearing the interrupt due to a leading 
edge on the control Ine is chosen. This routine not only 
clears the interrupt but also reads a timer and passes 
this value into memory and resets the tmer for the next 
leading edge. thereby the period of the signal on the 
control hne is turned into bytes, 

Tms mechnansm is expanded upon to give the data 
used by the processing program. The general fiow 
diagrams!2 for the data acquisition routines, *Acquis’ 
outline the approach in more detail. 

Acquis. interprets the signals into two groups of 
bytes by the process above and stores them to an inter- 
program buffer for use in a processing program, ‘Daps’ 
which finally gives the rig velocity and shaft encoders’ 
relative displacement!3, 

2.8 Problems and systematic errors 

There are two main hardware related errors generated 

by the interface as set up. One 1S concerned with the 


‘{see Figure 5, Interrupt servicing (adapte 


Colin; tnterfacing the BBC Microcomputer; 
1984) 

'2The flow diagrams may be found towards the 
document 


13gee Figure 6, Flow of infarmation 
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limit on the displacement measured to an eighteen degree 
arc and the other is due to the method by which the 


timers are read, 


2.8.1 Over 


_.wist 


The encoding dises are given a displacement of 
between four and nine degrees for the two reasons:- 

)..to allow the recording oF ‘twists: in” ‘the  direc- 

tion reverse to the direction of motion of the shaft, 

n)..to prevent the interrupt due to the start of 

the period from disc A coinciding with the interrupt 

for the next leading edge ont pene is zero or near 


zero twist. 


\f the twist on the shaft is more than the length of 
a black and white zone (or the arc remaining after the 
displacement mentioned above) then the leading edge of 
the signal from marking the end of ta would come after 


the next leading edge for the end of the period!4. 


The acquisition routines copes with iS problem by 
waiting to measure an interrupt of the correct nature. 
That is, if it has just processed and interrupt to start 
the period ~ and interval timers then it will not use 


another interrupt from that source until it recewes one 


14F igure hs Over twist of the encoder 
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that signals the end of the interval. The times 
transferred to the processing program refer to a double 


period. 


2.8.2 Timer Roll Over 


An error arises when reading the 16 bit wide timer 
eight bits at a time. The timer is not stopped when it is 
read and it takes a finite time to read even half of the 
clock, fo illustrate that this in certain circumstances 
iS a problem, consider a timer reading at the interrupt, 
ta) hexadecimal:- 

FoOe, 
the least significant byte, O02 is read into one of the 
crU registers, Ths takes three merosecongds, in which 
time the timer now reads:- 

FCFF 
as it counts down. The most significant byte, FC is now 
read, When the read value of the timer is sent to memory 
it would be recorded as:- 

Fco2 
precisely 256 (decimal) off the true value at the event. 
it is not possible for the acquisition program to guard 
against this. The time is sent as read to Daps, which 


corrects for these errors. 
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2.8.3 Variable Interrupt Latency 

A final error is due to the microcoding of the CPU 
used i microcomputers. there is a small variation in the 
times for the servicing due to the time taken to first 
detect the signal, the CPU only checks for the IRQ at the 
end of each machine code instruction the delay wn 
detection may be trom zero to two microseconds. this 
change lay delay is reflected wn the times read. 
2.9 Limitations of the system (BBC microcomputer) 
2.9.4 Frequency. of measurements 

(he processing time imts the frequency of inter- 
rupts the system can handle hence imits the frequency of. 
the signals. Ins hmitation constrains tne number of 
transitions marked on each of the encoding discs. The 
processing times for the service routines im Acquis when 
running on a BBC Master series microcomputer are given in 
the table below. 

the CRt cB2 

0.S 10 13.0 13.0 13.0 
entry 23.0 25.0 26.0 28.0 
service 09.5 20.5 40.5 23.0 
other < = 25.5 28.5 
exit o7.5 OTS OTS5 OF.'5 
TOTAL 53.0 66.0 S15 72.0 
where O.S. refers to the time for the operating system to 

pass the control to the service routines. 
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ventry’ is the time to decide on the particular service 
routine, 

*service’ is the time for clearing the condition and 
carrying out some specific action, 

’other’ in the case of the CB! and 2 routines is the time 
to send the received data to the man memory. 

‘exit’ the time taken to resume the interrupted = activity, 

T1, 2 and CBi, 2 refer to the source of the interrupt in 


Y the VIA. 


The maximum frequency of interrupts that the 


computer can cope with is given by the sum of the CBI, 


CBe service times and time for addition initial 
detection, that is: 
191.5 microseconds, 


The maximum frequency that may be dealt with is 
10,4 KHZ. 
Acquis places the data from the interface into the inter- 
program buffer when the CPU is not otherwise occupied by 
interrupts. The time required for complete transfer of 
data collected is 
321 microseconds, 
The time that must be allowed fur processing 1s 


435.25 microsecond 


This leads to the maximum frequency of interrupts for 
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complete transfer of 
2.3 KHz. 
A combined frequency 2@KHz is the maximum practicable to 


allow management of the buffers. 


2.9.2 Accuracy of time measurement 


With the BBC microcomputer the timers are clock at 


Mhz, There is a lhmt of one clock pulse on the accuracy 
of measurement of the times, Ths is a fixed absolute 
error, particular to the system, of 1 microsecond, As the 


periods and interval decrease in length with increasing 


speed the relative error also increases. 


Some digital velocimeters employ algorithms to 
reduce the relative error by the extension of the periods 
of measurement.’5 This approach is not open to the 
current approach as it reduces the frequency of 
measurements. The next section details improvements which 


reduce this absolute error. 


When dealing with a signal from a stmp made of 20 


black/white zones travelling at 1500 rpm, these residual 


'SBonnert, Richard, "Digital Tachometer with Fast 
Dynamic Response Implemented by a Microprocessor"; 
{EEE Trans Industry Applications 1A-19 (6); 
November/December 1983; pp 1052 - 1056 
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errors are: for the Master -series computers variation in 
latency 0-4 us 
clock 4 us 
accuracy 0.8% (0,16 degrees) 
It should be noted that the latency and clock perods are 
fixed so the accuracy will improve for a particular 
system when measuring at a lower velocity and 


correspondingly worsen at greater velocities. 


2.10 Improvements - Mark tl and beyond 


2.10.4 Transfer to a 32 bit RISC _ microcomputer 


With the introduction of a iow enat 32 obit) =6RISC 
based microcomputer from Acorn, the Archimedes, the 
opportunity was taken to transfer the developed software 
to this faster environment. A table above details the 
processing times for the routines of Acquis, On adapting 
certain routines the processing times within software 
have been reduced to one tenth that lsted, Ths left 
enough time for Daps to operate in real time even though 


it currently operates in interpreted Basic V, 


The ‘User Port’ suppled with this computer uses 
timers that are clocked at 2Mhz. This reduces’ the 
absolute error in the time measurements to 500 


nanoseconds. 
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To illustrate the limits of the hardware used: when 
dealing with a signal from a strip made of 20 black/white 
zones travelling at 1500 rpm, these residual errors are: 


for the Archimedes variation in 


latency 3.5-22.5 processor cycles 
2.375 us 
clock o.5 us 

= accuracy 0,3% (0,04 degrees) 


(Compare with the intervals for the Master given in the 


previous section) 


2.10.2 Mark u interface 

A large portion of Acquis is devoted to the 
manipulation of various registers within the VIA to 
generate the ee it iS proposed that most of these 
algorithms may successfully be transferred to hardwired 
logic, with the times presented to the CPU directly on 


its data bus. Circuit 2 iS a schematic diagram of part of 


the circuit. 


The advantages of this interface are the elimination 
of errors in reading due to variable interrupt - latency, 
the reduction in processing time per signal and the 
reduction of the absolute error in the measurements (by 


the use of a clock frequency »>10MHz), 


Using the Wiustration at the end of the previous 
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section, the error in measurements of twist will be 0.03% 


(0.004 degrees) 


2. Data Processing, Daps 


There are four routines that deal with the collection 
data, the correction for hardware faults and the final 
processing to velocity and shaft torsion, 


These procedures are named (in this report):- 


getbits {recovery of the dataj 

daterror {compensation for faults} 

process {calculation of the velocity and twist} 
-ratios {scaling factors for the coding strip} 


Except for ratios, every set of data has to be passed 
through these routines in turn to be = processed. They 
incorporate the algorithms described above, An 


interpretation in Pascal is given in Appendix us 


getbits 

This routine recovers the raw data stored by the 
acquisition program, acquis, nv the inter-program buffer 
convert them to integer times!®, Ihere are considerable 
differences between Master and Archimedes versions 
because of the differences ila the memory mapping and 


routines in Acquis. A Pascal version of this and its 


16, flow diagram for the Master version is given with 
the flow diagrams 
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subprocedures are given? although in practice machine 


code is used for speed. 


daterror 

There are a number of systematic errors) introduced 
into the collected data by the present interface that 
this routine 1s designed to detect and correct. These 


have been detailed im earher sections. 


- These errors are detected by this routine in the 


same way. Four point rolling averages of the periods, ty, 
and intervals, ta, are maintained by another routine 
process. Against the averages each new set of data is 
tested, 

if there is an error due to the over twist of the 


shaft during measurement then the period measured would 
be approximately twice that of immediately before. The 
period is nalved and time lapsed are corrected and 


processing continues. 


If over twist continues more than one period the 
interval would be much less than the average. If it falis 
to below one eighth of the previous value it is modified 


by adding the previous period. the following Pascal code 


‘Tappendix 1, getbits 
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ilustrates the principle elements of this routine. 


if sample<oldsample[1]/8 then 


sample:=samplet+oldperiod[1]; 

if (period>1,8*AVperiod) and (period<2.2x*AVperiod) 
then 
begin 


period: =period/2; 
runtime: =runtime+tperiod 
end; 

The ranges on the period are to bracket the differences 
in the times of the periods due to changes in speed. if 
the interval measured, the value of sample, is greater 
than one eighth of the previous measured it is taken as a 
true reading not subject to errors, By this method no 


correction is required when the shaft returns to the 


normal range of measurement. 


if period>AVperiod+150 then period: =period-256; 
if sample>AVsample+150 then sample: =sample-256; 

These. last two conditions of daterror detect the 
clock rollovers while reading, When either the period or 
interval (sample) are longer than the previous average by 
some large figure then they are reduced appropriately. 
The limit is set at 150 microseconds longer than the 
previous time because of the variation in speed and the 
coding disc strips but to allow the detection of the 


majority of these data faults. 
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process 


Ths updates 


rolling averages and recalculates them. 


velocity and code 


equations (4.1) and 


derived from the 


time represents 


runtime, and is used to end the iteration 


ratios 


This routine 


factors for the 


Unlike the above 


the main iteration 


their services. 


follows equations 


the points used 


disc displacement 
(4.2). For graphical 
accumulated periods 


the time within 


sets up the values 


individual zones 


mentioned routines 
loop of daps—= although 


To calculate the 


18see Appendix 1; 
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CHAPTER 3 
The Development Rig, Nature and 
Bet The Development Rig 
The test rig on which the monitoring technique 
been developed is described, it will be used further 
work into the nature of transients within induction motor 


drive systems 


The lst of components below 


exhaustive = but describe the = prin 


the test bed during the stage of 


ce The components of the 
Induction Motor 
Brook Motors Ltd 
3 phase induction motor 
Speed 
Coil Voltage 
Coil Current 
Frequency 
Output Shaft 
Number of rotor bars 
Number of stator tee 


Transmission Chain 
Bearings Fenners Ltd 
RHP 1020-20 Self al 


Shaft Mild Steel 

Pulleys Fenners Ltd 
SPZ type 
Dia. 
Bush 

Belt Fenners Ltd 


SPZ 1800 680 for 


19see Figure 8, Photograph of the 
schematic diagram of the same, 
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iS not intended 
ciple configuration 


the work 


development 


R10-A31910 
1440 rpm. 
400/440V 


1.25" dia. 
54 


th 36 


igning 
20mm dia. 


Cat No. 03120243 


160mm 
No, 20e1 


140/140 


rig and Figure 


has 


be 
of 


reported!?, 
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DC Generator 
Higgs Motors of Witton 
Model W1033-1 
Originally wired in compound mode 
Now with 
a separately excited field 
3.3 The Model 
The model consists of a system of first order 
ordinary differential equations for each part of the rig 
which are solved by the Runge-Kutta algorithm2° 


computer. €ach part of the model iS separately 


discussed, 


3.3.4 The induction motor@t 


The standard presentation of the electrical 
equations for a three phase induction machine involve a 
six by six matrix connecting the voltages and currents in 
the three phases for stator and rotor. The elements of 


the matrix are differential coefficients with time, 


fv } |r +#DL. OM DM OM DM OM ae 4 
vty Toh 1 re 4$BL) omt3— pmt4 opm? pm 16 vets 
ve{=}om@t pA 2 » @BL ome4 —pm25 m2 fw 2, 
v3} |om3! — pm32 pw 3 or #BL M35 pm 3e Aral 
v4l jom4! = pm42—pm43— oo 4 om4é ie 
v5} fom®1 = om®2 = pw&3 M84 ph re PBL | 6 83 
6 64 62 63 64 6 6 (3.1) 
20Thomas, Benku; "The Runge-Kutta Methods"; BYTE; 
4/1986; pp 191 - 210. 
@tuones, Charles V.; The Unified Theory of Electrical 


Machines; Butterworth, 1967 
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(3.1) may readily be deduced from the equivalent 


circuit of the induction motor but in this state it is 


cumbersome to solve. Jones®! offers transformations that 
will reduce the number of equations in the special case 
of balanced supply. it is hs model that is adapted for 


the prime mover. 


fi gph}=(Cal Lizpn). 
: (3.2) 


where {Ca} 1s the transformation matrix from the three 
phase currents, (iapn) to the resolved equivalent 
currents of a two phase machine, Clapnd- The 
transformation is square so the two phase current vector 
has six components of which two are zero one for the 
rotor system and the other for the = stator. By analysis 


it is found that the matrix is: 


c = BO 
C ro 
Oo B 
where B= (8743) 1/2 a 


ie) oO 
-(176) 172 (172) 142 Oo 
=(176) 172 ~(172) 1/2 0 


6323): 
With the transformation: 
(Zapn}= (Cal TiZ3ph} 
(3.4) 
and dropping the zero rows and columns one obtains: 
v2 PR aD {e} MOcos@ MDsin@ 
Rial Poy 1 R DL -MDsin@ MDcose | xt) i, 
[vr MDcos@ -Mbsind R +DL_ 0 enh 
v%@} IMDsin@ MDcos@ of © R 4+DL_ 
ye ar (3.5) 
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where {vj} are the voltages across the stator and rotor of 
the equivalent two phase machine and M is and equivalent 
mutual inductance between windings (or bars). 

The equations are still non-hnear as they contain 
terms that depend on the angular position of the rotor 
with respect to the stator,@ which depends on time. 

fo} = 4 is) ce) oO ' 

* o "4 ) ° 
ie) Q Sinwt coswt 
ce) ° coswt sinwt 
(3.6) 
where WwW is the rotor angular velocity and t the time. 


On applying [Gel: 


Cees R +L oD 30 i) MD A Te 
yO: OPN Sener n DD o ene, 
of “Mw MB Rieti) =i il) ' 
° Mo Mw LRo Reel Db iW 
Re ve d (3.7) 
The motor air gap torque is given by: 
T= (1/2) (13! (dae) CL My tng 
(3.8) 
Transforming {1} with [Ca] and (Cy) to give the d-q 
currents then (3.8) becomes: 


Rie 


From Jones®!, (d/d@){L,M) is 
-M |O oO oO sine 
DE to) io) sine 
° () (c) sine4 
sin@ sine sind Oo a 
sine sine4 sine® a 
singe sing sine4 o 
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sine 
sine 
sing 
oO 
oO 
oO 


(172) Figg)" LCg)! (Ca) §(d/d@) (LM) [Cal [Cg] Figg). 
(3. 


2 


4 


9) 


sing 
sine? 
sine® 
0 

oO 

Le) 
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where @5 = @+2i/3 
84 = O+40/3 
My 1S the mean value of the measured 


mutual inductance, 
(3.10) 
On substituting into (3.9): 
TMeM(igipDtigig@) 
where M is the value taken for the mutual inductance. 


(3.41) 


The currents for this equation are obtained from the d-q 


médel given wn equation (3,7) suitably manipulated to 


make the derivatives of current expheit, ready for 
numerical integration, 
Di L “LR, M@w ML Ww MR (1 
D 2D 2s 2 D 
Di Lov |+[-M@w -L Ro MR —ML_w |xfa 
a Ss a 
Di ~My ML w MR LR Lok ow 1 
q S$ R wove q 
Di ~Mv MR she Weel JL w Kc UR ' 
d d Ss 1 hire Nea: d 
(3. 42) 
3.3.2 the transmission chain 
3.3.2.1 Motor to Orive Shaft 


The drive Shaft iS assumed to be isotropic and 
homogeneous therefore is given a characteristic constant 
scalar torsional stiffness. The torque appled to the 


shaft causes an angular displacement of its ends, a 


twist, wn reaction, This is 
TsH = Koy (85p-85m)- 
Ss AA} 
The frictional and windage torques for the drive 
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shaft and the motor rotor are lumped and measured 
together. A cubic approximation angular velocity is used 
to model this drag. 

Try = Gy toaDAte3 (DO) %+e4 (00) 7 

(3.12) 

The equation governing the motion at the motor to 
shaft coupling, hence the velocity of the motor’s rotor 
and the displacement of this end of the shaft is:- 


D2@eq° (IMF Pet) Zl 4 
(3.13) 


3.3.2.2 Drive Shaft to Pulley 
Similar to the above the equation of motion of the 
pulley end of the shaft is given by 


D2agp=(IsH- Te) /la, 
(3,14) 


3.3.2.3 The | belt. transmission®2 
This requres a different approach from the = shaft 
transmission. A satisfactory has been gerived from a 
consideration of three aspects of the belt’s behaviour. 
i) The fractional speed loss between the pulleys; 
due to the normal sipping of the belt against the 
pulley surfaces while transmitting torque. 
(i) forque loss - the difference wn the torque 


applied at the driving pulley and that sensed at 


22n diagram of the belt and pulley is provided to help 
comprehension. Figure 10. 
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the following. this is taken as due to the 
frictional losses at the belts entry and exit pants 
to the pulleys. 

(ui) fhe global shp limt - The point at which the 
belt ceases to transmit torque in proportion to that 
apphed because it slides over the entire surface 


of one or other pulley. 


The model makes assumptions about the belt and 
pulley system that nave been demonstrated to be 
inaccurate by Cowburn®3, The model is not an accurate 
portrayal of the belts but gives a simulation of the 
transmission characteristics when embedded in the larger 
model. The short periods over which the model is_ solved 


allow the following approximations. 


Creep theory?4, where the belt 15 replaced by a thin 
elastic member passed over fiat pulleys, is the basis of 
the model. The effects that are considered the most 


import to model are detailed below. 


Speed loss 


The elastic belt stretches ia) proportion to the 


23cowburn, David; The mechanical performance of 
automotive V-belts; PhD. Thesis, University of 
Bradford; 1984 

24Gerbert, B. G.; "Force and slip behaviour in V-belts"; 


Acta Polytechnica Scandinavia Mechanical Engineering 
Series No. 67; 1972 
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tension set up within Kt. However, when transmitting 
torque there is a difference between the tensions in the 
two arms of the belt. This difference causes a torques on 
the pulleys ia) reaction, Hence the belt arms are nn 


different states of extension. 


AS the belt travels over the driving pulley it 
relaxes against the dynamic friction provided by the 
pulley surface, Ths friction supports the difference’ in 


tensions and the reaction torque. The change in extension 
causes the belt to have an exit speed lower than the 
pulley’s circumference, This process is reversed on the 


following pulley, 


In’ normal operation, the belt does not globally slip 
against the pulleys’ surfaces. There is at least one 
point on each pulley where the belt is held by = static 
friction and moves at the same speed as the surface, [hs 
point is at (or near) the first place of contact of belt 
and pulley in the direction of the belt velocity. The 
slack arm of belt is travelling faster than the following 
pulley and the taut arm slower than the driving pulley. 


The following pulley has the drive pulley’s belt exit 


speed and vice versa. So an output shaft rotates at a 
lower velocity than the input shaft. Itis is the speed 
loss, 
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Other effects operate on speed loss, the next most 
important being the radial comphance of the belt 
allowing it to ride out of the "V" groove used on this 
mg. 

From Gowburn®s, 

2s 


(D@gp-DOgp) /DOgp = QiFy-Fely 


q@ =. (17 © GYR). 


Torque _ toss 
There is some loss of the transmitted torque due to 
some. slipping at the entry and exit points where the 
beit is accelerated to the shape of the pulley, not all 
of the bending moment required at entry is recovered at. 
exit. ihe balance iS torque lost in friction, These and 
other effects are estimated to be, by Cowburne®3; 
Taree = 0.04 (140. 331n(105EI/R4K) I (F Et g)R 
for # £ 0.4, 
for wp » 0.4 substitute 0.4 for wp. 
Caer 
The torque loss acts to reduce the torque in the 
belt seen by the following pulley only and does not 


reduce the torque acting on the belt by the driving 


pulley, 


Global slip 
The belt slips with no pomt fixed relative to the 


pulley surface when the ratio of the tensions in each 
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arm exceeds the limit given by the capstan 


formula: 


Fi/Fs 2% Tpreexp (us) - 
(3.18) 


This may reduced to: 


(Ft-Fs)max = 2Tppotanntu/2es, 
(3.49) 


where (Fr.Fg)max 18 the greatest difference in 
tensions tenable between the two belt’s arms before 


giobal shppage. 


Model Equations 

From. equation (3.19), the afference in tensions may 
be written explicitly dependent on the instantaneous 
angular speeds of the driving and following, 


(Ft-Fs)max2 1/9) *(D8gp - DO¢gi/(ID@gp!), 
(3:20) 


Following the analysis: given by Jones, the torque 
supplied by the generator when loaded is given by 


Tg= DOG(NGSI FlatKla)s 
(3.24) 


with a correction for the remnant flux. 


This torque acts with the torque apphed from the 
pulley on the inertia for the generator rotor resulting 


in an angular acceleration, 
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D2aq = (Teg - To - Tral/!3, 
(3.22) 
where Tro is the frictional torque formulated similarly 


to equation (3.12), 


3.4 The mig model coefficients and parameters 

The detailed testing of the parameters for use 
within the model is described by the references for the 
model. The methods used follow Jones®! and CGowburn@s, 


Induction motor 
Electrical parameters 


Vie 420 volts 
R 4.5 ohms 

L 0.597 henrys 
M 0.565 henrys 


Number poles m3 


Mechanical parameters 
Vy 0.0421 Kgm@ 


for Ty 
cy  2,86x1073 
ca 3.78x1073 
cz -1.62x1073 


cg 3.47«1073 


Shaft and Belt transmission 


key 525 radiansN7!m-! 
R 0.09 m 

Kp 1.28%10°9 Nmo! 

lo 0.0382 kgm@ 


DC generator 
Electrical parameters 


Ye neglected 
kK O.14 
Gng 6.21 

rm 1.25 


No. poles 1 
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Mechanical parameters 


13 0.2239 kgm@ 
for Tro 

es Sarsxiors 

ce 1. 30x1074 


c7 2.56x1073 
cg -7.84x10-4 
ue 1.4 
3.5 The Model Equations in Summary 


THeMCigiptigiq)s 


where M is the value taken for the mutual inductance, 
(3.44) 
o Lov -L_R.  M@y~ = ML wR nee} 
tease tayo! almtake 2 Pale 
[Oi f= 4 {Usvp+[-M@@ -L RMR “ML ow [efit 
ee 2a 26 R. 2 Qa 
[OV | (LL o-M@)]}-Mv | {| ME@ MR -LR SC Th 
q jie = q 1 Ss 1R ‘ae q 
o ~My MRO -ML w -L Low -L_R ' 
ely Ear tet se (lie 
(3,12) 
TsH = KsH(@gp-9¢m) - 
; (3.44) 
Try = C4+G2D0+e3 (D8) “+4 (De) 3. 
(42) 
DF@sm= (Tu-TsH-Te1) 71 4- 
(3.43) 
D2agp=(Tsy- Tg) /la. 
(3.14) 
DT = 0,04 (140. 331M [105E1/R4K ,)}(FEtFs) Ry 
for » £ 0.4, 
for wv > 0.4 substitute 0.4 for p. 
pase =(1/qj«(0e - D@. 3/(}|be ‘ 
: t S max dp fd l api (3.20) 
(Ft-Fs)max = 2Tpretanh(ps2e). 
(3.149) 
D?0g = (Tpg - Te - Tra)/ la. 
(3.22) 
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3.6 Computer solution of the model 


The equations in section 3.3 are solved by a fourth order 


Runge-Kutta technique. 


Gy pthtp (ty Aye Kar tern Bry) 

Ga pthfp(ttH/2, 44464, 4/2, % 2404, 2/2, vee Kn tGq ,n/2) 
G3 ,ptHfp (trey %4+Ga, 4/2, *atG2, 2/2, teen Bn tGe2,n/2) 
G4, rat (t+, %1 +63, 44 KatG3, a4 vert ¥n*63,n) 
and 


Ap (at Tel) axpt [Gp 4 t2lGo, ptO3,p-}+G4,r]/6, 
where fp is the equation of the differential coefficient 


of Kp in terms of the state variables and time,25 


3.7 Ihe computer control of the rig 

Experiments will be carmed out on the states of 
twist that occur normally and due to ‘faults’ such as 
start-up, star-delta switching, fast bus transfer and 
various types of loading conditions, For repeatability 
and = flexibility the switching of the rig is controlled by 


thyristors Signalled by computer. 


3.7.4 Three Phase Supply Control 


The induction motor has two configurations of its 
stator windings to the three phases of the supply. These 


are ‘star’ and ‘delta’. The thyristor control is able to 


2Srully anotated flow diagrams are displayed in the appendices 


listings for the Acorn BBC Master Turbo are available 


BASIC IV and 1SO-Pascal, a tisting in BASIC Vv 
assembler for the Acorn Archimedes is also available 
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effect these two configurations and allow switching from 
one to the other26, The selection a particular 
configuration and switching transients reduces to that of 
the control of the sets of thyristors for each phase. 
fis cantrol iw provided by the firing circuts and the 


signals sent to the same from the computer interface®’, 


3.7.2 Generator toad 


Experiments mighhght the effects of loading 
variations on the angle of twist. The load on the 
Ganeceten is set by microcomputer by switching gate turn 
off thyristors (GTOs) that he wn parallel resistive 
circuits loading the generator armature28, The 
resistances are chosen to draw from @Q to 15.75 amperes 
from the generator in steps of O.2b5A, The control of the 


GTOs 183 Carried out by an interfaced firing eimourte9, 


3.7.3 Software Control 
Control of the motor switching and loading of the 
generator is written for the Archimedes only which has 


the capability cies interacting in real time while 


26see Figure 11, Star/Delta windings of an induction 
motor 

27see Circuits 2, The thyristor firing circuit; and 3 
The computer interface, 


28see Circuit 4, Load the DC generator 
29see Circuit 6, Interface to computer. 
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collecting data. All routines are written in RISC 
assembler code in the form of operating system primitives 
in a ‘relocatable module’ 39 to be used in higher level 


programs, 


39Acorn Computers, Archimedes Progrmmer’s Guide, 1988 


JtlI300988 49 


CNAA Registration report 


Chapter 4 


Outline of experimentation 


The experiments carmed out on the rig described 
above have two ams: the development of the software and 
mechanism of the twist detector to retrieve better grades 
of information and the andlysis of the effects of the 
variouS modes of operation of the mg in terms of the 
State of twist of the transmssion shaft. These two aims 
comeade "y the goal of continuous monitoring of this 


system from the twist signals, 


4a fhe start up transient 
This transient, which generates some of the severe 
torsional vibrations, iS used to configure the torsion 


meter and the supporting software, 


The intial results are compared to the expected 
result from the time domain model ot the rig. Any 
differences are carefully noted and the software 


inspected and improved, 


4.2 Other tests 

To discover parameters that distinguish between 
different transient and continuous modes of the system 
operation records will be taken of the state of twist for 


each of the cases. These will be analysed in both time 


o 
° 
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features. 


Comparisons will 


determine 
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Chapter 5 


The Results, Graphs from the Meter and Model 


2 Below are the results obtained from the torsion 


meter and the models of the mig. These are intended to be 
a clear example of the meter’s capabilities by a 
comparison of the values obtained with the theoretically 
predicted. 

These results have been collected by the Mark | 
interface on a BBC Master microcomputer. 

The final sections of the chapter explain the nature 
oF the errors that are corrected by and remain after the 


current data processing. 


5.1 Comparison between the experimental and theoretical3! 
A visual inspection of the closeness of fit 


between the model and the experimental results. 


A number of the electrical and mechanical 
aspects of the mg were monitored to gauge the accuracy 
of the lumped parameter model, These are presented as 
graphs 44 to a) and represent the unloaded and 


electromagnetically loaded start up transient. 


31 the graphs are drawn after the other figures 


a 
nm 
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5.2 Errors demonstrated by the time domain graphs 


6.2.4 Faithless recording of low speed transients 


The frequency resolution of the encoding discs is 
directly proportional to the speed. At full speed each 
encoding dise samples the velocity at approximately 
500Hz, This gives a maximum transient frequency of 2@50Hz 
recorded, However, at low speeds, 100rpm, the fastest 
recorded transient frequency becomes {7Hz,. The resonant 


frequency of the mg is about 2é@Hz, as is seen in the 
graphs, Compare Graphs 10 and 11. The first three cycles 
of twist are lost and the vibration is not clear until 


the encoding discs speed up. 


9.2.2 Zone5 of doubt ang sensitivity. 


Tnere are two zones at the top and bottom of the 
angular displacement graphs where the system would have 
to cope with the simultaneous arrival of two interrupts. 
Due to the interrupt latency the computer fais to record 
accurate data. The iwmaccuracy stems from the nature of 
the timers, In the period between the VIA issuing and the 
CPU acknowledging an IRQ the timer continues to count 
down. This problem will be discussed at length in Chapter 
7. it is this extension to the time recorded that causes 
the zones because the CPU would be delayed due to the 


immediately previous IRG. 
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6.2.3 Variable interrupt servicing 


jms oS a@ delay of between O and 4 microseconds on 
top of the time recorded due to the normal latency. The 
calculated error on the displacement record +/- 
0.08 degrees. 
— 

The width of the variation in the tal of the 
displacement Graph 10 is measured at 0.15 degrees, It is 


proposed that tho armiation iS due entirely to 


the variable delay im the interrupt service. 


5.3 Frequency domain graphs 

Only a limted amount of analysis of the frequency 
domains of the velocity and twist signals nas been 
earried out.3¢ the fast Fourver algorithm demonstrates 
the limitation of the system by the Nyqust frequency to 
approximately 250Hz. This constraint + due to the use of 


20 transitions per revolution for measurement. 


Some significant cumponents and there harmonics are 
apparent from the graph. These have been identified in 
frequency but Nave not been related sources in the rig or 


the processing. 


32560 graphs 12 and 1% 
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§.2.3: Variable interrupt servicing 


jus 18 a delay of between O and 4 meroseconds on 
top of the time recorded due to the normal latency. The 
calculated error on the displacement record +/- 
0.08 degrees, 

The width of the .variation in the tal of the 
displacement Graph 10 1S measured at O15 oR eee It is 


proposed that this. variation iS due entirely to 


the varrable delay ow the interrupt service. 


5.3 Frequency domain graphs 

Only a limted amount of analysis of the frequency 
domains of the velocity and twist signals nas been 
carried out.3¢ the fast Fourrer algorithm demonstrates 
the limitation of the system by the Nyqust frequency to 
approximately 250HzZ, This constraint i due to the use of 


20. 6transiticns per revolution for measurement. 


Some significant cumponents and there harmonics are 
apparent from the graph. These have been identified in 
frequency but have not been related sources in the mg oor 


the processing. 


32500 graphs 12 and 13 
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Chapter 6 


Concluding discussion 


Even at this stage when much expemmentation and 
analysis needs yet to be completed, signals, after 
elementary processing, are generated from the twist of 
the shaft that are sufficient sensitive and accurate for 
the purpose of deducing the state of the machine. There 
are strong indications that the frequency analysis will 
reveal a signature that gives indication of the torque 


signature and the mechanical transfer to the detector. 


AS mentioned these conclusions are as yet poorly 
founded and the remainder of the work will develop” the 
technique and explore the signals from the rig under 


different operating conditions. 


The work is expected to be continued to explore the 
application of this: technique to emore general and 


complex structures, for example a crankshaft, 


The basis of the thesis for the degree of Doctor of 


Philosophy will be developed from data collected from a 
crankshaft. There will be an exploration of the 
feasibility to determine the operating condition of an 
internal combustion engine by the signals from an 


encoder, as described here, attached to is crankshaft, 
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Appendix 4 


Sclected Routines in Puscal Ce 


procedure getbits; 
{recover data from the sideways ram for processing} 
constsraccess=66,; 
osword =&fff1; 
type word = 0. ,65535; 
fudgerec = packed record 
case boolean of 
false: (addr: “byte); 
true : (int :word ) 
end; 


var dummy: integer; 


function peek (address:word):byte; - « 
{emulates PEEK to inspect a memory location 
from Acornsoft guide to 1!SO-Pascal on BBC microcomputer} 
var fudge; fudgerec; 
begin 
fudge.int:=address; 
peek: =-fudge.addr® 
end{peek }; 


function samplereturn: integer; [turn data into an integer} 
var databyte, datavalue, jy: integer; 
begin 
datavalue: =0; 
for ji:=0 to 2 do 
begin 
databyte:=((-peek(dataty)-1) 
+ 256)MOD &100; 
(datavalue) *256 
+databyte; 


datavalue 


end; 
samplereturn: -datavalue; 
end{samplereturn); 


procedure bits;{get data from sideways ram} 

var dummy: integer; 
begin 
params[9):=readsample mod &100; 
params[10):=(readsample div &100) mod &100; 
dummy : =code1 (osword, sraccess, params) ; 
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number: =number +1; 
readsample:=-readsample+3; 
end; (bits) 


begin{body of getbits} 
bits; 
period; =samplereturn/1&b6; 
bits 
sample: =samplereturn/1&6; 
if (Number=2)and(samplerperiod) then 
begin 

period: ssample; 

bits; 

Sample: =samplereturn/te6 


end; 
end; [getbits} 


procedure ratios; 
var 1,g:integer; 


begin 
for i= to 1 
begin 
For hese. Co 20: 
begin 
Zone[i,j):=0; 
if ratiob then zoneli,JjJ=1; 
end; 
end; 
if ratiob then 
begin 
readsample:=12000; 
repeat 
begin 
for s:=f to 20 
begin 
getbits; 
daterror; 
process; 
zone[0,0):=zone[O,0)+period; 
zone[(O,1});=zone[(O,1)+period; 
zone([1,0):=zone[0O,1]+sample; 
zone(1,iJ:-zone[1,1] +sample; 
end; 
end 
until readsample>=36000; 
for 1:20 to 20 
begin 


zone(O,t)=zonel0,1)/zone[O,0)/5 
zone[1,1)=zone(1,1)/zone[1,0] /5 
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end; 
readsample:=0 
numberofsamples: =O . 
end; 
endjratios]; 


procedure daterror; 
{compensates for some hardware errors} 
begin 
if number>20 then 
begin 
if samplecsoldsample(ij/8 then 
sample: -sampletroldperiod{1}; 


if (period71,g*AVperiod) and (period<2,2*AVperiod) 


then 
begin 
period: =period/2; 
runtime: =runtime 
end; 


eriod-256 
ample-256; 


if period>AVper:od+150 then period 
if Sample>AVsample+150 then sample 
end; {daterror} 


procedure process; 
{turn the integer times into velocity and twist} 
var yrinteger; 
begin 
if ratiob else 
begin 
Sample:=samplexzone[\, 
(numberofsamples/d+offset)mod2o0 +1) 
period: =period«zone[o, 
(numberofsamples/2toffset)mod20 +1) 


end, 

for ji=3 downto 1 do 

begin 
oldsamples[jt1):=-oldsamples([yjJ; 
AVsample :sAVsamplet+oldsamples(Jj); 
oldperrods[(yt!) oldperiods[ jj; 
AVperiod AVperiod+oldperiods[y]) 

end; 


oldsamples[1j:=sample; 
oldperiods[(1)i=period; 
AVsample:=(AVsample+sample) /4; 
AVperiod: = (AVperiod+period) /4; 
if AVprocB then 
begin 

period: =AVperiod, 
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sample; =AVsample 


end; 
velocit O/period/stripe; 
runtime untime+period; 


lwist:=Sample/period 
end; (process} 
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Figure 1, Code disc location 
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Figure 2, Code disc construction 
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Figure 5, Interrupt servicing 
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Figure 7, Over twist of encoder 
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Figure 9, Schematic of rig 
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Figure 11, Star/delta connection of motor 
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General flow diagram of acquis 
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Circuit 1 - Opto-switch to TTL interface 
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Circuit 2 - Thyristor firing circuit 
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Circuit 4 - GTO controlled loading of the generator 
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Graph 2: Velocity - no timer rollover 


The routines for the detection and correction are 
applied to the data for graph 1 
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Graph 3: Velocity - as previous and no over-twist 


The routines for the detection and correction of over 
twist are applied to the data 
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Graph 6: Zones of doubt 


The areas of the graphs where interrupt latency 
interferes with the accuracy of the measurements 
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Graph 9:Angular Displacement - as previous and no over twist 


| 
| 
| 


1.60 1.80 


1.40 


41243 


Wh 


0.80 1.00 


Mon,O1 Feb 1988.19 


0.60 


0.40 


} 


| 
| 


Displece 
18.00 
14.00 
12.00 
10.00 | 
8.00 
6.00 
4.00 
2.00 
0.00 


jt1300988 86 


| 


Dieplacens 


CNAA Registration report 


Graph 10:Angular Displacement - as previous and disc 
compensation 
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Graph 12:Frequency analysis - velocity, of the uniform 
velocity section 


(0) 250 


FFT af 192, samples 


Jti300988 89 


CNAA Registration report 


Graph 13:Freq. analysis - twist, uniform velocity section 
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Graph 14:Verification of model - acceleration, motor 
model 


Loaded electromagnetically 


ts — + ae ———— _ 


2.00 Thlme/seconds 


0.80 1.00 1.20 1.40 1.60 1.80 


0.60 


2.164 
1.5E3 
9. 0€2 
-3.0E2 
-9.0EZ 
-1.5E3 
=2,1£ 
-2.7E3 


3.0E2 


Jjt1300988 os 


FIWATANARF INSTRUMENTS CORP. 


CNAA Registration report 


Graph 15:Verification of model 
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Graph 17: Verification of model - velocity, motor, 
experiment 
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Graph 417: Verification 
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Grapn 19:Verification of model - current, motor, 
experiment 
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Graph 21:Verification of model 
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